Relativistic recoil corrections to the atomic energy levels 

V. M. Shabaev 

Department of Physics, St. Petersburg State University, 

Oulianovskaya Street 1, Petrodvorets, St. Petersburg 198504, 

Russia 

(—^ ' The quantum electrodynamic theory of the nuclear recoil effect in atoms to all orders in aZ 

f^ , and to first order in m/M is considered. The complete a.Z-dependence formulas for the relativistic 

rvj ' recoil corrections to the atomic energy levels are derived in a simple way. The results of numerical 

calculations of the recoil effect to all orders in aZ are presented for hydrogenlike and lithiumlike 
atoms. These results are compared with analytical results obtained to lowest orders in aZ. It is 
ry^ ' shown that even for hydrogen the numerical calculations to all orders in aZ provide most precise 

theoretical predictions for the relativistic recoil correction of first order in m/M. 
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I. INTRODUCTION 



In the non-relativistic quantum mechanics the nuclear recoil effect for a hydrogenlike atom is easily taken into 
account by using the reduced mass /i = mM / {m + AI) instead of the electron mass m [M is the nuclear mass). It 
means that to account for the nuclear recoil effect to first order in m/M we must simply replace the binding energy 
E hy E{1- m/M). 

^ , Let us consider now a relativistic hydrogenlike atom. In the infinite nucleus mass approximation a hydrogenlike 

c/3 ■ atom is described by the Dirac equation {h = c= 1) 

C/2 . (-ia • V + /9m + yc(x))V'(x) = eVW, (1.1) 

>>■ 
1-^ ' where Vc is the Coulomb potential of the nucleus. For the point-nucleus case, analytical solution of this equation 
yields the well known formula for the energy of a bound state: 

^^, (1.2) 
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^j^ ' where n is the principal quantum number and j is the total angular momentum of the electron. The main problem 

f~^ • we will discuss in this paper is the following: what is the recoil correction to this formula? 

f^ It is known that to the lowest order in aZ the relativistic recoil correction to the energy levels can be derived 

^^ from the Breit equation. Such a derivation was made by Breit and Brown in 1948 |]1| (see also g|). They found 

c/3 that the relativistic recoil correction to the lowest order in aZ consists of two terms. The first term reduces the fine 

O , structure splitting by the factor (1 — m/Af ). The second term does not affect the fine structure splitting and is equal 

c/3 ' to —{aZYm^/l^Mn^). Calculations of the recoil effect to higher orders in aZ demand using QED beyond the Breit 

^^' approximation. In quantum electrodynamics a two-body system is generally treated by the Bethe-Salpeter method 
Q or by one of versions of the quasipotential method proposed first by Logunov and Tavkhelidze W. In Ref. ^J 



(see also B), using the Bethe-Salpeter equation, Salpeter calculated the recoil correction of order {aZym? /M to the 

energy levels of a hydrogenlike atom. This correction gives a contribution of 359 kHz to the 2s - 2pi/2 splitting in 
L^ ■ hydrogen. The current uncertainties of the Lamb and isotopic shift measurements are much smaller than this value 
r\ , (see, e.g., 0) and, therefore, calculations of the recoil corrections of higher orders in aZ are required. In addition, for 
j^ the last decade a great progress was made in high precision measurements of the Lamb shifts in high-Z few-electron 

ions ||q-[l(]|]. In these systems, the parameter aZ is not small and, therefore, calculations of the relativistic recoil 

corrections to all orders in aZ are needed. 

II. RELATIVISTIC FORMULA FOR THE RECOIL CORRECTION 

First attempts to derive formulas for the relativistic recoil corrections to all orders in aZ were undertaken in 0,n2|. 



As a result of these attempts, only a part of the desired expressions was found in |12| (sec Ref. |l3| for details). The 



complete aZ-dependence formula for the relativistic recoil effect in the case of a hydrogenlike atom was derived in [Il4| 



The derivation of [Q was based on using a quasipotential equation in which the heavy particle is put on the mass 
shell 1 15 1^ . According to ||lj] , the relativistic recoil correction to the energy of a state a is the sum of a lower-order 
term AEi^ and a higher-order term A£'h: 

AE = AEi^ + AEu, (2.1) 



A£;h = / duj 



xG(c^ + £,)(DM + lPl^)|a). (2.3) 

V uj + lO / 

Here \a) is the unperturbed state of the Dirac electron in the Coulomb field Vc{r) — —aZ/r, p ~ — iV is the 
momentum operator, G(w) = [lo — H{1 — «0)]~^ is the relativistic Coulomb-Green function, iJ = a • p + /3m -I- Vc , 
ai (l — 1,2,3) are the Dirac matrices, Sa is the unperturbed Dirac-Coulomb energy, 

Drn{uj) = -ATraZaiDijr,{uj) , (2.4) 

Dimiijj) is the transverse part of the photon propagator in the Coulomb gauge. In the coordinate representation it is 

n I r.\ 1 f cxp(ii|cj|r) (exp(ii|cj|r) -1) \ 

Dik[ijJ,v) ^ ---< Oife + ViVfc 5 \. (2.5) 

The scalar product is implicit in the equation ( |2.3| ). In Refs. |17| , [l8| , the formulas ( |2.l|) -( py3|) were rederived by other 
methods and in |17| it was noticed that AE can be written in the following compact form: 

AE^^-^J dco{a\[p^-Diu;)]G{Lu + ea)[p--Diu;)]\a). (2.6) 



However, the representation ( |2.lD -( p.3D is more convenient for practical calculations. 

The term A_El can easily be calculated by using the virial relations for the Dirac equation ||l^,Q . Such a calculation 
gives |l|l 



m^ — e^ 



AEi^ = ^ . (2.7) 

2M ^ ' 

This simple formula contains all the recoil corrections within the (aZ)^m^/Af approximation. The term AE'h taken 
to the lowest order in aZ gives the Salpeter correction H. Evaluation of this term to all orders in aZ will be discussed 
below. 

The complete aZ-dependence formulas for the nuclear recoil corrections in high Z few-electron atoms were derived 
in Ref. |£l|. As it follows from these formulas, within the {aZy-m^/M approximation the nuclear recoil corrections 
can be obtained by averaging the operator 

s,s' ^ 

with the Dirac wave functions. This operator can also be used for relativistic calculations of t he nu clear recoil effect 



in neutral atoms. An independent derivation of this operator was done in |22[| . The operator (2.8) was employed in 
p3[ to calculate the {aZ)^m'^/M corrections to the energy levels of two- and three-electron multicharged ions. 

III. SIMPLE APPROACH TO THE RECOIL EFFECT IN ATOMS 

As was shown in p3| , to include the relativistic recoil corrections in calculations of the energy levels, we must add 
to the standard Hamiltonian of the electron-positron field interacting with the quantized electromagnetic field and 
with the Coulomb field of the nucleus Vc, taken in the Coulomb gauge, the following term 



Hm = ^ydx^^(x)(-zV,)^(x)y dyV'^(y)(-jVy)^(y) 



p7 r p^7^ 

■— / dxV.nx)(-*V.)^(x)A(0) -I- ^^A2(0) . (3.1) 



This operator acts only on the electron-positron and electromagnetic field variables. The normal ordered form of Hm 
taken in the interaction representation must be added to the interaction Hamiltonian. It gives additional elements 
to the Feynman rules for the Green function. In the Furry picture, in addition to the standard Feynman rules in the 
energy representation (see |gj,^), the following verteces and lines appear (we assume that the Coulomb gauge is 
used) 

1. Coulomb contribution. 

An additional line (" Coulomb-recoil" line) appears to be 

^% r du. 

2-K M J— oo 

X y 

This line joins two vertices each of which corresponds to 



W3 



-27rj7°(5(cji - UJ2- ^z) f dupk 

where p = — jVx and k — 1,2, 3. 
2. One-transverse-photon contribution. 

An additional vertex on an electron line appears to be 



W3 



W2 




-2TTij°6{uJi -UJ2 -W3)^/dxpfc , 
The transverse photon line attached to this vertex (at the point x) is 

X y 

At the point y this line is to be attached to an usual vertex in which we have — 27rie7°a/27r5(a'i — cj2 ~ ^a) / dy, 
where a/ {I — 1, 2, 3) are the usual Dirac matrices. 

3. Two-transverse-photon contribution. 

An additional line ("two-transverse-photon-recoil" line) appears to be 



2tt M J— oo 



/-^ ^'^^«('^; x)i:)ifc(w, y) ■ 



This line joins usual vertices (see the previous item). 

Let as apply this formalism to the case of a single level a in a one-electron atom. To find the Coulomb nuclear 
recoil correction we have to calculate the contribution of the diagram shown in Fig. 1. A simple calculation of this 
diagram yields (see Ref. jl^] for details) 

The one-transverse-photon nuclear recoil correction corresponds to the diagrams shown in Fig. 2. One easily obtains 



_ AnaZ i f°° ^^ J (a|pj|n)(n|afcAfe(ea - ^)|a) 



a\akDtk{£a-uj)\n){n\p^\a)_] 



The two-transverse-photon nuclear recoil correction is defined by the diagram shown in Fig. 3. We find 

(47raZ)2 i 
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{a\aiDu{£a ~ Lo)\n){n\akDik{£a ~ ^)\a) ,„ .-. 

"" c.-e„(l-*0) ■ ^^-^^ 

The sum of the contributions (^)-(^) is 

1 i f°° 
^^=T7^/ duj{a\{p, + ATTaZaiDu{uj)) 
M Ztt 7_^ 

xG{ijj + ea){pi + 47raZa,„Am(w))|a) ■ (3.5) 



This exactly coincides with formula (p^ 

Consider now a high-Z two-electron atom. For simplicity, we will assume that the unperturbed wave function is a 
one-determinant function 

u(xi,X2) = — =^(-l)^f/;pQ(xi)V'P6(x2). (3.6) 

The nuclear recoil correction is the sum of the one-electron and two-electron contributions. The one-electron contri- 



bution is the sum of the expressions (3.5) for the a and b states. The two-electron contributions are defined by the 



diagrams shown in Figs. 4-6. A simple calculation of these diagrams yields 

A^('"*) = ^ Y^{-lf{Pa\p^ + AT:aZaiDu{ePa - £a)|a) 
p 
X {Pb\p, + inaZarnDrmiepb - Sb)\b) . (3.7) 

The formula (3/7) was first derived by the quasipotential method in pi[ . 



IV. NUMERICAL RESULTS 
A. Hydrogenlike atoms 



According to equations (2.1)-(2.3) the recoil correction is the sum of the low-order and higher-order terms. The 



low-order term AEj^ is given by equation (|2.7| ). The higher order term AEn was calculated to all orders in aZ in 



5kE7|. The results of these calculations expressed in terms of the function P{aZ) defined as 



M -Kv? ^ ' 
are presented in Table 1. To the lowest order in aZ the function P{aZ) is given by Salpeter's expressions: 



>(i«) 



P. 



{2s) 



{aZ) = --log {aZ) - - 2.984129 

2 8 

(aZ) = --log (aZ) - - 2.811769 ^ 
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Comparing the function P{aZ) from Table 1 with the lowest order contributions ( [4.2[ )-(4.4) shows that for high Z the 
complete aZ-dependence results differ considerably from Salpeter's ones. 

In the case of hydrogen, the difference AP = P — Ps amounts to -0.01616(3), -0.01617(5), and 0.00772 for the Is, 
2s, and 2pi/2 states, respectively. Table 2 displays the relativistic recoil corrections, beyond the Salpeter ones, to the 
hydrogen energy levels. These values include also the corresponding correction from the low-order term ( p.7[ ) which 
is calculated by 



A'£;[^^) 
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(4.5) 
(4.6) 



The results of Refs. |25|j27|| which are exact in aZ are compared with the related corrections obtained to the lowest 
order in aZ. In [ p8| , ^9[ it was found that the (ctZ)^ log (aZ)m? / M corrections cancel each other. The [aZ^'ir? jM 
correction was derived in [Q for s-states and in [gO| for p-states. The (a.Z'f\o'^{aZ)rr? jM correction was recently 
evaluated in Refs. |31 33]. The uncertainty of the calculation based on the expansion in olZ is defined by uncalculated 
terms of order (aZymr jM and is expected to be about 1 kHz for the Is state. It follows that the results of the 
complete aZ-dependence calculations are in a good agreement with the results obtained to lowest orders in aZ but 
are of much higher accuracy. 

As it follows from Ref. [O], the formulas (2.1)- ( pT^ ) will incorporate partially the nuclear size corrections to the 
recoil effect if Vc(f) is taken to be the potential of an extended nucleus. In particular, this replacement allows one to 
account for the nuclear size corrections to the Coulomb part of the recoil effect. In Ref. p^], where the calculations of 
the recoil effect for extended nuclei were performed, it was found that, in the case of hydrogen, the leading relativistic 
nuclear size correction to the Coulomb low-order part is comparable with the total value of the [aZ^^vr? jM correction 
but is cancelled by the nuclear size correction to the Coulomb higher-order part. 

One of the main goals of the calculations of Refs. 12^,^270 was to evaluate the nuclear recoil correction for highly 
charged ions. In the case of the ground state of hydrogenlike uranium these calculations yield -0.51 eV for the point 
nucleus case |25 and -0.46 eV for the extended nucleus case |Q. This correction is big enough to be included in 
the current theoretical prediction for the Is Lamb shift in hydrogenlike uranium p4| but is small compared with the 
present experimental uncertainty which amounts to 13 eV |10| . However, a much higher precision was obtained in 
experiments with heavy lithiumlike ions 0,0]. In this connection in Refs. |25 23] the nuclear recoil corrections for 
lithiumlike ions were calculated as well. 



B. Lithiumlike ions 



In lithiumlike ions, in addition to the one-electron contributions, we must evaluate the two-electron contributions. 
In the case of one electron over the (Is)^ shell the total two-electron contribution to the zeroth order in XjZ is given 
by the expression 



A£;'"' = -^ ^ (a]p-D(ea-£„)]n)(n]p-D(ea-e„)]a), 



(4.7) 



where D is defined by equation (p]4). Calculation of this term causes no problem ]p3,pq]. For the 2pi/2 a-nd 2^3/2 
states, the results of this calculation expressed in terms of the function QiaZ) defined by 



2^777^ 

Ai?'"* = -33— (aZ)2Q(aZ) 



(4.8) 



are presented in Table 3. For the s-states the two-electron contribution is equal zero. To the lowest orders in aZ the 
function Q{aZ) is given by [^ 

Q['''^^\aZ) = 1 + {aZ)'{-% + log^) , (4.9) 

Q';;'^^^\aZ) = 1 + {aZr{-'^^ + ilogf) . (4.10) 



The expressions (4.9)-( 4.10 ) serve as a good approximation for the Q{aZ) function even for very high Z . 

For low Z , in addition to the corrections considered here, the Coulomb interelectronic interaction effect on the non- 
relativistic nuclear recoil correction must be taken into account. It contributes on the level of order {l/Z){aZYm? /M . 

To date, the highest precision in experiments with heavy ions was obtained for the 2p^/2 — 2s transition in lithiumlike 
bismuth ||]. The transition energy measured in this experiment amounts to (2788.14±0.04) eV. In g the energy of the 
2pi/2 ~ 2s transition in lithiumlike uranium was measured to be (280.59 ± 0.10) eV. In both cases the recoil correction 
amounts to -0.07 eV and, therefore, is comparable with the experimental uncertainty. At present, the uncertainty 
of the theoretical predictions for these transition energies is defined by uncalculated contributions of second order in 
a (see Refs. [^,Q). When calculations of these contributions are completed, it will be possible to probe the recoil 
effect in high-Z few-electron systems. This will provide a unique possibility for testing the quantum electrodynamics 
in the region of strong coupling {aZ ^ 1) beyond the external field approximation since in calculations of all other 
QED corrections in heavy ions the nucleus is considered only as a stationary source of the classical electromagnetic 
field. 

V. CONCLUSION 

In this paper the relativistic theory of the recoil effect in atoms is considered. It is shown that the complete aZ- 
dependence calculation of the recoil correction provides the highest precision even in the case of hydrogen. The recoil 
corrections to the energy levels of highly charged ions contribute on the level of the present experimental accuracy. It 
provides good perspectives for testing the quantum electrodynamics in the region of strong coupling {aZ ^ 1) beyond 
the external field approximation. 
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FIG. 1. Coulomb nuclear recoil diagram. 
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FIG. 2. One-transverse-photon nuclear recoil diagrams. 



FIG. 3. Two-transverse-photon nuclear recoil diagram. 



FIG. 4. Two-electron Coulomb nuclear recoil diagram. 



FIG. 5. Two-electron one-transverse-photon nuclear recoil diagrams. 



FIG. 6. Two-electron two-transverse-photon nuclear recoil diagram. 



TABLE I. The results of the numerical calculation of the function P{aZ) for low-lying states of hydrogenlike atoms. 



Is 



2s 



2pi/2 



2p3/2 



1 


5.42990(3) 


5 


4.3033(4) 


10 


3.7950(1) 


20 


3.2940(1) 


30 


3.0437(1) 


40 


2.9268(1) 


50 


2.9137(1) 


60 


3.0061(2) 


70 


3.2334(4) 


80 


3.672(1) 


90 


4.519(8) 


100 


6.4(1) 



6.15483(5) 

5.0335(2) 

4.5383(1) 

4.0825 

3.9037 

3.8900 

4.0228(1) 

4.3248(2) 

4.8656(5) 

5.807(2) 

7.557(9) 

11.4(2) 



-0.30112 

-0.2692 

-0.2277 

-0.1393 

-0.0421 

0.0685 

0.2000 

0.3655 

0.5894 

0.9214(2) 

1.481(1) 

2.63(2) 



-0.3013(4) 

-0.2724(1) 

-0.2379 

-0.1726 

-0.1107 

-0.0517 

0.0050 

0.0597 

0.1125 

0.1638 

0.2138 

0.2625 



TABLE IL The values of the relativistic recoil correction to hydrogen energy levels beyond the Salpeter contribution, in 
kHz. The values given in the second and third rows include the {aZ)^m^/M contribution and all the contributions of higher 
orders in aZ. In the last row the sum of the {aZY'in^ /M and (aZY log^ {a Z)m^ /M contributions is given. 



State 



Is 



2s 



2p 



1/2 



To all orders in aZ, Ref. [25] 
To all orders in aZ, Ref. [27] 
{aZfm^/M, Refs. [18,30] 
laZy\og\aZ)m^/M, Refs. [31,32] 
The sum of the low-order terms 



-7.1(9) 
-7.16(1) 
-7.4 
-0.4 

-7.8 



-0.73(6) 

-0.737(3) 

-0.77 

-0.05 

-0.82 



0.59 

0.587 

0.58 



TABLE in. The results of the numerical calculation of the function Q{aZ) for low- lying states of lithiumlike ions. 



(l^)"2pi/2 



(l^)"2p3/2 



10 
20 
30 
40 
50 
60 
70 
80 
90 
100 



0.99741 
0.98959 
0.97645 
0.95776 
0.93313 
0.90195 
0.86320 
0.81529 
0.75570 
0.68041 



0.99810 
0.99239 
0.98281 
0.96926 
0.95165 
0.92988 
0.90390 
0.87362 
0.83896 
0.79951 
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